Gravothermal Catastrophe with a Cosmological Constant 
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We investigate the effect of a cosmological constant on the gravothermal catastrophe in the New- 
tonian limit. A negative cosmological constant acts as a thermodynamic 'destabilizer'. The Antonov 
radius gets smaller and the instability occurs, not only for negative but also for positive energy val- 
ues. A positive cosmological constant acts as a 'stabilizer' of the system, which, in this case, exhibits 
a novel 'reentrant behaviour'. In addition to the Antonov radius we find a second critical radius, 
where an 'inverse Antonov transition' occurs; a series of local entropy maxima is restored. 



I. INTRODUCTION 



Antonov's gravothermal instability is an important ef- 
fect in gravitational thermodynamics [THl] . It has served 
as the prime paradigm for an extensive research into the 
statistical mechanics of systems with long range interac- 
tions in different fields of physics [H, @ , while phase tran- 
sitions of self-gravitating systems is an active research 
field 0. One may ask if and how this phenomenon de- 
pends on the vacuum background. A positive cosmolog- 
ical constant is nowadays one of the main candidates for 
dark energy [f| |j| , while negative cosmological constant 
has attracted a lot of attention due to the anti-dc Sit- 
ter/ Conformal Field Theory correspondence (AdS/CFT 
[lOJ ) . In addition, the stability properties of anti-dc Sit- 
ter spacetime have become recently a subject of topi- 
cal interest [ill, EH ■ We present here how gravothermal 
catastrophe is affected by a cosmological constant, i.e. by 
a de Sitter or anti-de Sitter vacuum. For convenience we 
call the non-relativistic limit of de Sitter [l3[ and anti-de 
Sitter spaces (usually called Newton- Hooke spaces [3]), 
just dS or AdS, respectively. 

The system under study in the original formulation 
[l|, [|J is a self-gravitating gas in the Newtonian limit 
bound by a spherical shell with insulating and perfectly 
reflecting walls with fixed energy (microcanonical ensem- 
ble) and fixed number of point particles (stars). An 
equilibrium state corresponding to an entropy maximum 
is called an 'isothermal sphere'. Antonov proved that 
for such a system in the mean field approximation there 
is no global entropy maximum. However, there ex- 
ist local entropy maxima (i.e. metastable states) for 
E-R > — 0.335GM 2 . This means that for positive energy 
E, there are isothermal spheres for every radius R, but 
for negative energy local entropy extrema exist only for 
radii smaller than a critical value Ra, we call Antonov 
radius. For R > Ra and fixed negative energy there does 
not exist any equilibrium state. 

For a fixed negative energy and for radii smaller than 
the Antonov radius R < Ra, the equilibrium state may 
be stable (local entropy maxima) or unstable (saddle 
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points) depending on the value of the ratio of core den- 
sity po = p(0) versus the edge ratio pr = p(R). There is 
a critical value {po/ pr)ct = 709 at which an instability 
sets in, i.e. for po/pn > 709. A remarkable feature of 
gravity is that this instability in the microcanonical en- 
semble sets in when the specific heat goes from negative 
to positive values and not the other way around (non- 
equivalence of ensembles). 

We find that gravothermal catastrophe (described in 
the last two paragraphs) depends crucially on the vacuum 
background. AdS space destabilizes the system. Com- 
pared to the flat case, the instability sets in at smaller 
radius, at higher central density and occurs not only for 
negative but also for positive energies. In dS space the 
phenomenon is drastically altered. A series of equilib- 
rium solutions is restored for large radii, indicating a type 
of 'reentrant' behaviour. A metastable homogeneous so- 
lution that suffers a transition to Antonov instability is 
found, together with metastable states that do not suf- 
fer any transition to instability. In addition, centrally 
diluted and periodically condensed solutions are allowed. 



II. ANALYSIS 



Consider a self-gravitating gas of N point particles, 
each of mass rh = 1, inside a spherical shell with insu- 
lating and perfectly reflecting walls in the presence of a 
cosmological constant [l5|. We work in the mean field 
approximation and use the one body distribution func- 
tion f{r,p). This is viable as long as the correlations 
between the particles are not significant [!, [ljl [ljj • The 
form of the distribution function at the equilibrium is 
found by the extrcmization of the Boltzmann entropy 
S = —k J / log fd 3 rd 3 p under the constraints of con- 
stant energy E and number of particles N (i.e. constant 
mass M = N ■ fh). It is easy to generalize for AdS and 
dS vacua, Antonov's proof that a global entropy max- 
imum does not exist. Only local entropy extrema do 
exist. In Ref. [l], [|| is proved for the fiat case that 
only spherical distributions maximize the entropy. In 
this work we consider only spherical distributions fl8j |. 
We find that at an equilibrium the distribution function 
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is: f(r,v) = (/3/27r) 3 / 2 e -^ 2 / 2 p(r) , where 

p(r)=poe-W-«°» (1) 

is the density of matter and 4>{r) is the potential that 
obeys the equation [l9| 

V 2 = AnGp - 8ttGpa (2) 

where the density p\ is related to the cosmological con- 
stant by the equation pa = Ac 2 /87rG. The potential can 
therefore be written as cj> = <f>N + 4>h with: 

(j) N = -G J |-r— — d r , (pA= ^ PisX (3) 

Using the spherical symmetry and the dimensionlcss vari- 
ables y = /3(cp — (f>{0)), x = ry/iTrGpof] and A = 2p^/po 
equation ([2]) becomes: 

which wc call the 'Emden-A' equation. This equation 
indicates that for every p\ > 0, there exists a radius 
Rh for which a series of homogeneous solutions, with 
p = 2/9a = const., <j)(r) = <f>'(r) = 0, exists. This radius 
depends only on p\ with Rh = (independent of 

energy). The homogeneous solution resembles the Ein- 
stein's static universe in the non-relativistic limit. 

In general, we want to solve the Emdcn-A equation 
with initial conditions y(0) — y'(0) = for various 
isothermal spheres, characterized by the radius R, keep- 
ing M constant. Let call z = R^/AnGpoft the value of 
x at R. The cosmological constant introduces a mass 
scale M\ = ^nR 3 p\ and we define the dimensionless 
mass to = M/2M\. Using the dimensionlcss tempera- 
ture (3 — GM/3 1 R, to can be written as 

to = 3/3/Az 2 (5) 

Dimensionless temperature can be calculated at z, by 
integrating the Emden-A equation: 

[3(z) = zy'(z) + l -\z 2 (6) 

In order to keep to fixed for various z we see from (J5]) 
that A has to be different at each z. We developed an 
algorithm [2(| that calculates the right A for each z and 
solves Emden-A with respect to z, keeping to constant. 
Writing m as to = J — Hr? it is clear that solving for 
various fixed m can be interpreted as solving for various 
PA and/or R for a fixed M. 

Using the virial [H [H theorem: 2K + U N - 2U A = 
3PV we derive the following expression for the non- 
dimensional energy at z, Q(z) = RE/GM 2 : 

^, ^ z 2 e-y^ 3 A r 4 _ v(x) , 

Q(z) = =- = =— / a; 4 e y(x) dx (7) 

y ' (3 2 2/3 2f3 2 zJo 1 ; 
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FIG. 1: In dS, the dimensionless energy Q = RE/GM 2 ver- 
sus log(po/pi?) for fixed M, R, p\. Every point of any curve 
corresponds to a distinct isothermal sphere, so that a curve 
represents a series of solutions. Unlike flat case, this plot 
cannot be interpreted with varying R and fixed E, because at 
each point, m is held constant. There exist many other curves 
not shown at each case for more negative energies. Only the 
upper series correspond to p(r) monotonically changing. Dis- 
tinct series corresponding to more negative energies have more 
extrema of p(r). At points B\, B2, -B3, -B5 and Ba an instabil- 
ity sets in, while curves A3B3 and A5B5 are already unstable 
solutions. 



Using ([7]) Q can be calculated numerically. For dS it is 
drawn in Figure [TJ 

In order to conclude on the type of stability we use 
Poincare's linear series of equilibria theorem 0, Oil [2l| . 
It states that as S(E) varies from one equilibrium to an- 
other by infinitesimal changing the parameter E (with E 
being held fixed in order to calculate the extremum of S) , 
the type of the extremum does change if this linear series 
of equilibria meets another series or turns back through 
previous values of E. The later is our case. Therefore 
at an extremum of E, the stability changes. Comple- 
mentary to this criterion we use the second variation of 
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entropy (following closely Padmanabhan Q): 



d 2 S 



R r R 



8M(r2)K(ri,r 2 )5M(ri)dr 1 dr2 
where SM is a local mass perturbation and: 



(8) 



K(n,r 2 )= - 
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\ Anpr 2 dr\ ) 



(9) 



The sign of 5 2 S is determined by the sign of the eigen- 
values £ of the eigenvalue problem: 



K(r,r 1 )Ft(r 1 )dr 1 =ZFe(r) 



(10) 



with Fg(0) = F^(R) = 0. We developed an algorithm 
that can numerically determine eigenvalues and eigen- 
states (the perturbations) of equation (JTUJ). If, for an 
equilibrium, equation (|10[) has solutions with one or more 
positive eigenvalues £, then this equilibrium is unstable, 
because for this mode S 2 S > 0. 

The instability of an equilibrium can be proved by find- 
ing just one unstable mode. However, the stability is 
much harder to be proved. A strategy is to prove stabil- 
ity for a limiting case and claim that the stability does 
not change across various E until E reaches an extremum 
(Poincare). We verified numerically this claim for the 
specific solutionspresented here. In the flat case, stabil- 
ity is proved P, Q for R — > 0. In this limit the presence 
of a cosmological constant is irrelevant. Thus, the same 
argument holds for AdS and for dS only when R < Rh 
for the curve A\B\ of Figure 1(a) If R > Rh this argu- 
ment does not hold and we prove stability for the curve 
A4B4 of Figure 1(b) as follows: the limit (3 — > does 
exist for these series of equilibria 1 since it corresponds to 
E -> 00. We see that S 2 S < for fi -> 0, when 6 2 S is 
written in the form: 



S 2 S = - 




Gf3_ 



{SM)' 



dr 



(11) 



8npr 2 

Note that in equations © , (jlip the cosmological constant 
enters implicitly through the derivative of the potential. 

For the homogeneous solution at R = Rh we find 
that although E has not an extremum, there exists a 
critical point S 2 S = determined analytically, with 
Zcr = 4.493 (3 cr = 6.73, where the stable series of equi- 
libria become unstable. It is evident that such a point 
should exist, since /3 — > leads to a stable state (equation 
(JTT|)), while (3 — > 00 gives the unstable static state. 



1 Intuitively, one can argue that for ft — > the kinetic energy 
dominates, therefore the system behaves as an ideal gas, hence 
it is thermodynamically stable. 
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FIG. 2: The critical energy and critical ratio of the centre 
to edge density for fixed radius R and fixed mass M vs the 
cosmological constant. (a)For E < E cr (unshaded region) 
there is no equilibrium state, while in the shaded region there 
exist equilibria (stable or unstable, depending on the ratio 
po/pB.)- (b)For po/pR < (po/pfl)cr the equilibrium is stable, 
while for pa/ PR > (po/pr)ct it is unstable. 



III. RESULTS 

For a fixed radius R < Ra, where metastable states 
with negative energy do exist, we find that the critical 
energy E cr and the critical ratio (po/pr)ct ar e decreas- 
ing with increasing cosmological constant p\ as one can 
see in Figure [2] E cr is the minimum energy for which 
an isothermal sphere exists. For AdS the instability 
at (po/pR)cr sets in not only for negative, but also for 
positive energies, because the potential energy of mat- 
ter coupled to AdS is positive. In Figure [3] can be seen 
the critical radius R cr versus the cosmological constant 
Pa- As AdS becomes stronger, the region of no equilib- 
rium (R > Ra) is getting bigger. All these facts further 
support the recent investigations of the newly discovered 
non-linear instability of AdS spacetime [ill, El ■ 

The dc Sitter case is more complicated, though very in- 
teresting features arise. The critical Antonov radius Ra 




FIG. 3: The critical radius R cr for a fixed negative energy 
E and mass M versus the cosmological constant pA- For a 
given large enough and positive p\ there exist two critical 
radii Ra and Ria- In the white unshaded area there exist no 
equilibrium states. The line R — Rh is the radius for each 
pA where a homogeneous p — 2p\ — const, solution exists. 
In the region I there exist equilibria for which p(r) is mono- 
tonically decreasing that are metastable or unstable states, 
depending on the ratio po/pR. In the region II there exist 
metastable states for which p(r) is monotonically increasing 
that do not suffer a transition to Antonov instability. The 
small gray shaded region contains equilibria with p(r) not 
monotonic with one or more maxima and p(r) > po. This 
region stops at p A lm = p/4. 



increases with increasing pa defining the blue 'Antonov 
line' in Figure [3] There also appears a second line, where 
an 'inverse Antonov transition' occurs (black line in Fig- 
ure [3j) defined by a second critical radius Ria- This 
radius is decreasing with increasing pa- The two crit- 
ical radii merge at an extremal value p™ ax . This is 
a typical 'reentrant 1 behaviour. Reentrant phenomena 
of this type are common in statistical systems, when- 
ever competing interactions are present [22h25| . For a 
fixed pa, for R < Ra metastable states do exist and 
for Ra < R < Ria no thermodynamic equilibria exist. 
However, for R > Ria the system 'reenters' to a phase, 
where metastable states are possible. dS acts as a sta- 
bilizer for the system, since as dS becomes stronger the 
region of no equilibrium diminishes. For pa > p™ ax there 
are metastable states for every R and fixed E. 

One may intuitively understand the existence of a sec- 
ond critical radius Ria as follows. It is this radius beyond 
which the cosmological force that is exerted on the outer 
parts is so powerful that can hold enough mass at the 
outer regions, so that the remaining matter in the inner 
regions shall not collapse. 

As already mentioned in the analysis, the homogeneous 
solution is stable for temperature T > T cr , while for 
T < T cr is unstable, with T cr = 6 ■ Apart from the 
homogeneous solution, for R = Rh there exist infinitely 



FIG. 4: Schwartzschild-dS space has two horizons, a black 
hole horizon Rbh and a cosmological horizon Rc- In this 
figure, the horizon radius Rh is measured in units of the 
Schwartzschild radius Rs = 2 ^* f and is plotted w.r.t. the 
cosmological constant for a fixed mass. 



many other series of equilibria with various configura- 
tions p(r). These with p(r) < po arc unstable, while 
these with p(r) > p are stable. 

For R < Rh, in region I of Figure O there exist 
metastable states which correspond to the line A\B\ of 
For these, p(r) is monotonically decreas- 



1(a) 



Figure 

ing and can become unstable depending on the value of 
Pol PR likewise flat case. In the small upper gray shaded 
region of Figure [3J there exist equilibria with p(r) not 
monotonic, with one maximum and p(r) > po, V?*. These 
correspond to the curve A^Bi of Figure 1(a) Numeri- 



cal calculations indicate that this curve corresponds to 
stable solutions. At the maximum of energy B2 an insta- 
bility sets in. There exist a lot of other series with more 
maxima of p(r) that are not drawn in Figure [T] All these 
solutions correspond to isothermal spheres diluted at the 
center and with concentric condensations away from the 
center. 

At the region R > Rh (region II in Figure [3]) there ex- 
ist metastable states with p(r) monotonically increasing 
(curve A4B4 in Figure 1(b)), that do not suffer a tran- 



sition to instability. These isothermal spheres are some- 
what hollow with the mass being mostly concentrated 
at the edge. In addition there exist all other solutions 



shown in Figure 1(b) 



The small upper gray shaded region of Figure [3J con- 
tains the solutions A2B2 of Figure 1(a) along with other 



similar solutions for more negative energy that cannot 
be seen in Figure 1(a) All these, are finite in number 



for a fixed p\. Their multitude is diminishing as p\ is 
getting smaller. They exist only for pa > p™ m with 
purvm _ pj^ where p is the mean density. Therefore, the 
whole small gray region of Figure |3J stops at p™ m and 
its lower boundary does not cross the orthogonal axis. 



The value 



Pa 



can be calculated as follows. The 
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density contrast of series A2B2 behaves as po/pn — s- 
for diminishing cosmological constant. The p™ m is this 
value, for a fixed radius and mass, for which the cosmo- 
logical force is so strong that marginally can hold the 
entire mass at the edge, i.e. GM/2R 2 < ^-p h R => 
PA > 3M/16tt.R 3 so that p^ ln = 3M/16irR 3 = p/A. This 
result is supported by the numerical calculations, as well. 

Let us discuss the validity of the Newtonian limit with 
a cosmological constant (see Appendix, as well). In or- 
der to meet for A the weak field condition, required by 
the Newtonian approximation, the following condition is 
imposed when taking the Newtonian limit in Einstein's 
equations with A. 



1 



-Ai? 2 < 1 



(12) 



It is straightforward to verify that this condition leads to 
the condition 



„ IRs 

m ^> 

2 R 



(13) 



where m = 3AI/8irp\R 3 is defined in equation ([5]) and 
Rs = 2GM/c 2 is the Schwartzshild radius of the sys- 
tem. Condition (fT3"]) defines the validity of Newtonian 
approximation. The value of m that corresponds to the 
extremal p™ ax of Figure [3] is m ~ 1.57, satisfying per- 
fectly condition (TIB"]) , along with all other points of R cr 
of Figure [31 since they have m > 0.9. Another issue that 
might be raised is whether in general, A is negligible for 
Newtonian systems. That is not true. From Poisson with 
A equation ([2|) we see that the cosmological constant is 
negligible for a Newtonian system, if 



1.9 R 
P > Pa & xAi?- • — < 1 

6 Hs 



(14) 



Apparently, equation (fT"2|) docs not necessarily imply 
equation (|14j) . It may very well happen that condi- 
tion ([12]) holds with AR 2 ■ (R/Rs) ~ 1 for systems with 
R/R s > 1. 

An example from the physical world, when condition 
(TT2l holds but not condition (TTJ} is galaxy clusters and 
supcrclustcrs. There is evidence that the cosmological 
constant plays a role to the evolution of the galaxy clus- 
ters and superclusters [HI, [26| . In addition gravothermal 
catastrophe on its own, could play a role, as well, since 
many present a core-halo structure and others have a 
supcrmassive black hole in their centres. For some typi- 
cal values [H, [27} of M ~ 10 15 A/ Q , R ~ 2 - 5Mpc, E ~ 
— 10 57 J of regular galaxy clusters, their mean density p is 
of the same order of magnitude as the cosmological con- 
stant ~ 10~ 26 kgr/m 3 . Thus 2p^jp, that is the variable 
of the horizontal axis of Figure [3J is of order unity. For 
the same typical values, the quantity 8ttG 3 M 5 pa/3\E\ 3 , 
i.e. the variable of the horizontal axis of Figure [3j is of 
order unity, as well. This is evidence that our analysis 
could be relevant to the evolution of galaxy clusters. Wc 
believe that the effect of gravothermal catastrophe to the 



evolution of galaxy clusters on its own, and secondly the 
possible impact of the cosmological constant on it, de- 
serve further investigation. 

Last but not least, let us stress out that the reentrant 
behaviour of Figure [3] resembles the two horizons of the 
Schwartzschild-dS black hole in General Relativity. The 
radii of these horizons arc plotted in FigurcUwith respect 
to the cosmological constant. The similarity between Fig- 
ure [3] and Figure [4] is too striking to be considered a co- 
incidence! It seems that the reentrant phenomenon we 
have discovered is the closest Newtonian analogue to the 
two horizons of Schwartzschild-dS space. However, there 
is a big difference. The stable region of the Newtonian 
case corresponds to the unstable region of the relativistic 
case. This opposite sense, with no cosmological constant 
present, is explained in Ref. (28j . 
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Appendix 

We analytically derive from General Relativity the Em- 
den with A equation and discuss the validity of the New- 
tonian approximation. 

The Emden equation is well known to be the New- 
tonian limit of the Tolman-Oppenheimer-Volkov (TOV) 
equation for dust (non-relativistic) particles. We will 
show that similarly the Emden-A equation is the Newto- 
nian limit of the TOV equation with A. 

The relativistic TOV equation in the presence of A is: 



dp 
Jr 



~f= - (P/C 2 +P) 



GM(r) 



2GM(r) 



rcr- 



with 



dM(r) 



dr 



: Anpr 2 



4ttG 
—z-p-r 

c z 

-1 



-Ac 2 r 



(A.l) 



(A.2) 



For a perfect fluid the equation of state is p = wpc 2 with 
w = S-tt. where rh is the rest mass of one particle. For 
non-rclativistic particles it is w — ¥ 0, i.e. 

kT <C rhc 2 

Substituting equation p — wpc 2 into the TOV-A equation 
(|A. II) and taking the dust (non-relativistic) particles limit 
kT rhc 2 , we get 



kT dp _ (GM{i 



rh dr 



\kc 2 r 



(A.3) 
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Note that the term (2GM(r)/r)kT/rhc 2 is negligible so 
that in general relativity with A = just the limit w — > 
gives the Emden equation [29[, provided there is no sin- 
gularity. 

In order to obtain the Newtonian limit one has to as- 
sume that the gravitational field related to the cosmolog- 
ical constant is weak. We see from the above equation 
(|A.3I) that this accounts for: 

-AR 2 < 1 , Newtonian limit (A. 4) 

3 

since r < R. In this limit and substituting p\ = 
Ac 2 /8nG, j3 = 1/kT equation (|A.3|) becomes 

d[— insrlogp) 1 
r 2 \ gm *M = _ i_ 8 ^ GpAr 3 (A 5) 

dr S 



which after differentiation w.r.t r and substituting equa- 
tion (|A.2[) gives finally: 



1 d 

r 2 dr 



dr 



AnGp - 8irGp\ 



(A.6) 



with p = Ae P™^. This is the Emden- A equation. It is 
evident that the cosmological constant is negligible if 



1 



R 



p > p A ^ -Air ■ — < 1 , negligible A (A.7) 
3 Hs 

where Rs = 2GM/c 2 is the Schwartzscild radius. The 
equation (|A.4[) for the Newtonian limit does not nec- 
essarily imply equation (|A.7[) for negligible A since for 
Newtonian systems it normally holds that R/Rs 3> 1. 
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